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The other symbols in the text are defined as when they are 
introduced. 



ABSTRACT 


Load flow sLudy Is very important, for Lhe purpose of 

planning, operation and control of power system both for off-line 
and on-line applications. EXjring past four decades, extensive 
research efforts have been ejcpended in the area of load flow 
studies. Amongst various load flow methods available, Newton 
based load flow methods as well as fast decoupled load flow 
techniques have been widely used by utilities. Newton -Raphson 

CNR2> methods in polar as well as in rectangular co-ordinates are 

well known for offering reliable convergence to load flow 
solutions whereas fast decoupled load flow method in polar 
co-ordinates CFDLPZ), commercially developed, is known for less 
memory and CPU time requirements as compared to the two NR 

methods. Though, the convergence behaviour of these three methods 
are well established for well behaved systems, it was felt 
important to observe their performances for ill-conditioned 
situations as well. Hence, the motivation behind the work 
presented in this thesis, is to observe critically the relative 
F>erformance of the two NR methods and the FE4..P method in solving 
load flow problems for both well behaved as well as 
ill-conditioned cases in the systems. 

Moreover, literature survey shows that very little attention 
have been paid to the development of FES-F method in rectangular 
co-ordinates. An attempt has, therefore, been made in the present 



thesis to develop some new algorithms tor fast decoupled load Tlow 
method in rectangular co-ordinates CFDLRD and to compare their 
relative peformances with the existing methods. 

In the present work* computer programs have been develop>ed to 
perform load flow studies by existing Newton-Raphson & fast 
decoupled load flow methods In polar as well as in rectangular 
co-ordinates and by the new FEM_.R algorithms. A comparative study 
of all these methods with reference to computation time» number of 
iterations and memory requirement for both normal and 
ill-conditioned situations in the system have been presented. 
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CHAPTER - I 


INTRODUCTION 


1 . 1 Gener al : - 


Th* «1 metrical jsowor sys'Lctm ie continuously expanding in size 
and growing in complexity all over the world. Extensive system 
studies on computers^ both at off-line planning stage as well as 
at on-line operational stage» are therefore, a must. A power 
system engineer is always concerned with adopting newer and better 
mathematical models for solving system problems more efficiently 
and effectively. Load flow study is one such subproblem, which is 
so vital for almost all power system studies. 

A large number of research efforts have been expended over 
the past four decades in the area of power flow Cload flow!) 
studies and yet there continues a constant and persistent search 
for better and still better solution algorithms for load flow 
study in terms of speed, memory and reliability aspects. This is 
due to the fact that the load flow study is of fundamental nature 

V 

and any small improvement would be of considerable Importance to 
the electric utilities and pow&r companies. 

The main objective of the work presented in this thesis is to 
carry out in-depth study on some existing load flow methods and to 
explore the possibility of developing new load flow algorithms. 

Load flow study is a steady state solution of the flower 
system problems which provides information about network power 
flows and bus voltages. Load flow studies on a digital computer 
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&I I y Involves the solution of e set of non-linear algebraic 
equations. 

During the past forty years, probably almost all the relevant 
known methods of numerical analysis for solving a set of 
non-linear algebraic equations have been applied in developing 
load flow algorithms. In addition, many methods based on some 
spiecial features of the load flow equations have also been 
suggested. The JiK^st recurrent question arising in the load flow 
field is - which is the best method to choose for a given 
application? The answer is rarely easy. The relative properties 
and performance of different load flow methods can be influenced 
substantially by the type and the size of problem to be solved, by 
the computing facilities available and the precise details of 
iirqpilementation. The various desirable features of load flow 
methods can be : 

- High speed, specially for real-time applications 

Low and compact storage, specially for small computers 

- Reliability of convergence both for well-behaved and 
ill -condi ti oned systems 

- Simplicity and ease of programming. 

None of the methods available satisfy all the above criteria 
completely. Fortunately, not all the desirable features of load 
flow methods are essential in all the situations. A particular 
combination of the various features is needed in a given situation 
and hence, relentness effort continues to develop new load flow 
methods which would be nK»re suitable for certain types of systems 
than the existing methods. 



1.2 Li t.<»rat.ur«r Survey and Mot,i vaLlon : 


With th* enormous amount, of literature available in the area 
or load flow study ► it is not possible to present an exhaustive 
survey of all of them. Moreover ► an excellent review paper by 
Sftott [16] provides a detailed account of the load flow methods 
developed till 1974. Further, some of the common load flow 
algorithms are well documented in a number of books [2S, 26, 27, 

255, 293. Hence, the survey attempted here mainly confines to 

brief and critical review of some of the important works and 
especially those relevent to the investigation carried out in this 
thesi s. 

Some of the earliest methods developed for load flow solution 
were based on successive displacement iteration schemes. Amongst 
them, the two nodal Iterative methods Cl, 2] are based on bus 
admittance matrix CY 3 formulation also popularly known as <3auss 

BUS 

- iterative and Gauss - Seidel iterative methods. The attractive 
features of these methods are their simplicity and minimal storage 

- requirements. They are easily programmable and do not require a 
close estimate for initial bus voltages as is the case with most 
of the Newton based methods. However, these methods were slow in 
convergence and thus not very well suited to handle large systems. 

The slow convergence of iterative method is mainly 

because of highly sparse nature of, bus admittance matrix causing 
loose mathematical coupling between the buses. With this in view, 
load flow methods [3,5,83 based on bus matrix CZ 3 formulation 

tttISS 

were developed and found to take lesser number of iterations 
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and v/ara mora r«rllabl« than tha [Y 3 iterative models. However* 

■us 

the presence of P-V buses adversely affects their convergence 
behaviour. Further, because of increased computational time and 
storage requirements for computing and storing non-sparse * 

the method did not gain widespread acceptance with the utilities. 

An attempt was made by Freris and Sasson COl to combine the 
desirable features of EY 3 and [2 3 iterative methods by 

BUS BUS 

having first few iterations by algorithm and switching over 

tc algorithm for subsequent iterations. However, this 

method in general suffers from the same drawbacks as of pure 
[2 3 model . 

BUS 

A number of algorithms based on the powerful Newton-Raphson 
CNR!) method having quadratic convergence propierty have been 
proposed. The NR method essentially Involves the linearization 
of the set of load flow equations at any iteration point using 
only the first order term in Taylor’s series expansion. Several 
versions of NR method have been proposed and the px^ipular ones 
are based on power mismatches. The power mismatch versions can be 
formulated using either the polar E63 or rectangular co-ordinates 
C43 for voltages. Some other versions are based on current 
mismatches C73. The power mismatch version in polar co-ordinates 
has advantage over rectangular version in that it requires lesser 
memory space as number of variables are reduced. . The current 
mismatch version works well only in lightly loaded systems. The 
power mismatch version in polar co-ordinates developed by Tinney 
and Hart [63 utilizes the sparsity exploitation and optimal 
ordering of buses and has been widely used by the utilities. This 
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met,hc>d usually converges t,o a fairly accurate solution in 4 to 5 
iterations for most of the systems independent of their size and 
kind. Since, the convergence of the NR method is known to 
depend on good initial estimates of the variables otherwise there 
will be a large amount of truncation error, one or two iterations 
of Gauss -Seidel method or some other methods Cl 23 may be used 
before applying NR method for load flow solution. Constant 
Jacobian version of the NR method Cl 1,1 S3 and some other 
modifications [133 have also been reported. 

Some attempts have also been made to use minimization methods 
[10,203 for solving load flow problems. A minimization of sum of 
squares of current or pjower mismatches as objective function is 
considered whose global optimal coincides with the load flow 
solution. The minimization techniques used are gradient method 
[103, and linear programming method [203. As concluded by Barr as 
et al . [203, these methods are not competitive with NR algorithm 

for base case load flow application but are limited to contingency 
analysis. 

Though N R load flow methods proved to be highly reliable, 
the storage requirement is considerable higher as compared to 
<3*uss -Seidel iterative methods. Decoupled versions of NR 
method [14.153 in polar form were attempted in view of reducing 
memory and computational time. In all these methods, the weak 
coupling between the active power -vol tage magnitude and reactive 
power -vol tage angles were e^quiloited and two sets of linear 
equations, one to determine the corrections for voltage angles and 
the other to obtain the corrections for voltage magnitudes, were 
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solved in each iteration. Amongst various decoupled versions 
fast decoupled load flow CFDI-F^ method developed by Sftott ani 
Alsac 1153 is the most popular with utilities. The method is ii 
F»ola.r co-ordinates and utilizes certain Justifiable network 
assumptions apart from the P-V and Q-6 decoupling to obtalr 
constant Jacobian submatrices. The final model of FE4-F is 
developed employing certain additional assumptions in computing 
elements of the submatrices CB*3 and CB"3. The FEO^* model 
developed is quite fast and reliable for a wide spectrum of 
systems. Wu C173 and Rao et al . [1Q3 have tried to develop some 

empirical relations to study the convergence behaviour of FDLF 
method and have found it to exhibit pxxDr convergence in certain 
cases of ill-conditioned systems. Nanda et al . [S3] have examined 

the relative weight age and importance of various additional 
assumptions made by Stott and Alsac in their FDLF model on its 
convergence behaviour. Van Amerongen [243 suggested a different 
combination of considering additional assumptions in FDLF model 
which is more suitable for solving all types of systems. 

From the literature survey presented » it is observed that 
Newton -Raphson [4>6] and fast decoupled [153 methods are the most 
propular ones for load flow solution. Newton Raphson CNRD 
methods in polar as well as rcn^tangular co-ordinates are well 
known for offering reliable convergence to load flow solutions 
whereas fast decoupled load ,flow method in polar co-ordinates 
CFDLPJ is known for less memory and C. P. U. time requirements as 
compared to the two NR methods. Though the convergence 

behaviour of ifjhese three methods are well established for well 
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behaved sys'tems, it was felt important to observe their behaviour 
for ill-conditioned situations in the systems. Hence> the 

motivation behind the works presented in this thesis was te 
observe critically the relative performances of the two NR methods 
and the FCA-P method in solving load flow problems from the view 
point of memory requirements^ number of iterations and C. P. U. time 
for both well behaved and i 1 1 -eondi ti ons situations in the 
systems. 

Moreover > literature survey also shows that very little 
attention has been piaid to the development of FESLF method in 
rectangular co-ordinates. Only two attempts seem to have been 
made by Arumugam [213 and Singh et al . [223 in this direction, 

which incidently results into similar formulations. An attempt 
has therefore also been made in the present thesis in the 
development of some new algorithms for fast decoupled load flow 
method in rectangular co-ordinates CFDLRD and compare their 
relative p>erformances with the existing methods. 

1.3 Outline of the Thesis 

The chapter -wise summary of the works, covered out in this 
thesis is given as follows 

Chapter -2 is devoted to the theoretical and numerical aspects 
of the NR method in rectangular co-ordinates CNRRD and also the NR 
method in polar co-ordinates CNRP3. Both of these methods are 
tested on IEEE 14-Bus and 30-Bus test systems. A brief account of 
each method and their relative merits are also discussed in this 


chapter . 



Ch^p-tArr-S d«fiils wi'th Tasit. d«couplftd load Tlow nwrlhods In 

polar as well as in rectangular co-ordinates. The two new 
algorithms for FOLK are also presented in this chapter. Potential 
of all these methods have been tested on IEEE 14-Bus and 30-6us 
test systems both for well-behaved and ill-conditioned cases. The 
performances of various versions of fast decoupled load flow 
methods have been compared. 

Finally* the Chapter -4 eoncludfes the specific findings in 
this thcrsls along with future scope of the work. 



CHAPTER - II 


NEWTON BASED LOAD FLOW METHODS 


2.1 In-Lroduetlon 


The objective of load flow study is to determine the voltage 
magnitude and its angle at each bus, real and reactive power flow 
in each line and line losses of the p>ower system network for 
spyecified bus or terminal conditions. The mathematical model of 
the p>ower system network for the purpose of load flow studies are 
a set of non-linear algebraic equations which can be solved by 
Iterative techniques such as 6au&s method, Gauss-Seidel method or 
Newton’s method etc. 

The Newton-Raphson’s method is a powerful method of solving 
non-linear algebraic equations. It works faster, the number of 
iterations required to obtain solution is nearly independent of 
problem size and kind and is sure to converge in most cases as 
compared to the Gauss-Seidel CGS> method. In this chapter, 
Newton-Raphson’s CNR!) method in rectangular co-ordinates and also 
in polar co-ordinates for solving load flow problem are presented. 
The potential of both these Newton based methods for solving load 
flow problem have been critically examined for well behaved as 
well as ill-conditioned situations in the systems. 

t 

^ Load Flow Technique ; 

The load flow studies are conducted for the purpose of 
planning Cviz. short, medium or long range planning!), operation 
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and con'troil . For the purpose of load flow studies* we assume that 
the three phase power system network is balanced and we also 
neglect mutual coupling between elements. Variables associated 
with each bus of the p>ower system include four quantities viz. 
voltage magnitude V* its phase angle <5, real power P and 
reactive F»ower Q, Out of these four quantities* only two are 
generally specified at a bus and depending up»on which two are 
specified* we have following three categories of buses : 

aZ> Load or a P-0 Bus : 


For this type of bus* we know apriori load values CP^ ^ and 
O^ ^ 2> and specify the generations CP 
specify the bus injections P^^ and O^. Solution of the load flow 
equati ons will r ender and angl e 6^ . A 1 oad bus whi ch * due to 

its lack of generating equipment* is characterized by zero P^^ and 
also falls in this category. Here P and O &re specified. 
Sometimes V limits are also spyecified. 

bZ) P-V or a voltage controlled bus : 

For this type of bus* we know apriori P^ ^ ^ and specif yV^, 

and P^. . We thus specify P. . Solutions of the load flow 
equations render Cand hence and angle <5^. Such type of 

buses are also called a voltage controlled bus because its voltage 
can be controlled. 'Ibus at P-V buses* real power P and voltage 
magnitude V are specified alongwith Cl“limits to hold voltage 


Gi 


and 


^i' 


In fact we thus 


within limits. 



c5 Slack or cwl ng buss 


1 ] 


This is ■tha r«rf<rr«rnc6 bus whara t-he voll-age magni'tude V 
and phase angle S are specified. Voltage phase angle of such bus 
is normally sp^ecified as zero, thus it acts as reference to other 
bus voltage vectors. Usually a generator bus with ma>d.muin real 
power caf>ability is selected as a slack bus. The variables of 
interest at this bus are the real and reactive powers. Generators 
at this bus is also required to supply additional real and 
reactive piowers to meet the transmission losses which are unknown 
until the final load flow solution is obtained. 

2. 3 Newton -Raphson Method in Rectangular Coordinates : - 


Assuming balanced S-phase conditions as mentioned above, the 
transmission system network can be represented by its p»ositive 
sequence network. The nodal admittance matrix CY 1 which can be 
found by inspection, neglecting mutual coupling between elements. 

can be expressed for a *n' bus system as follows : 


tl 3 . = CY 3 [V 3 . 

BUS nxl BUS nxn bus nxl 


C2. ID 


The above equation can be written in the following form for 
the pth bus 


- I 


CY V D 

p pq q 

q=l 


C2.2D 


n. 


where, p = 1 ,2 



Thiss equa-tiori simply slat^fs thaLl curr^pnt at &ny bus i 

the 4silgebraie sum of all the currents entering or leaving the bus 


TTie power at any bus *p* CP^ - calculated using th 

relation : 


n 


P - JQ = V** I = V** V CY V 5 

p p pp p / ■ pq q 


C2. 3D 


q=l 


In rectangular coordinates, we express 


V = e + jf 
P P P 


Where e^ and are the real and imaginary components of the 
bus voltage at bus *p* 


Y = G - JB 

pq pq pq 


C2. 4D 


Where G and B are conductance and susceptance part of the 

pq pq 

pqth element of matrix respectively. Thus eqn. C2. 3D can 

be rewritten as. 


n 

-•p - J'V = '"P - 

q=l 


- JB D Ce + Jf D] C2. 5D 

pq pq p q 


From eqn. C2. 5D, we determine the net injected real and 
reactive powers at bus p. 


P 

P 


Real tv* 
P 


.ri— 'I 


Y .VI 

pq q 


n 


= I 


» G 

q pq 


+ f B D+f Cf G 

q pq p q pq 


q=l 


e B 

q pq 


D3 


C2. GD 



And 


% 


I magi nary 



I 


Y V 3 
p<? q 


n 

Z tr C® G + f B 5 - ® Cf G -® B 53 
pqpq qpq pqpqqpq 

q=i 


ca. 75 


Equat/ions Ca.G5 and Ca.75 at. ®ach bus tom t.h® 1 oad/''pov/c>r- 
flow equations which involve four variables P»Q, e and f at each 
bus. Thus* for an- bus system* total number of unknowns are 4n 
whereas the number of available equations are Bn only. Solution 
of these equations is obtained by specifying the values of two 
unknowns Cout of four 5 at each bus and thus reducing the total 
number of unknowns to an only C equal to the number of equations to 
t)e &olved5. Variables* which are usually specified at various 
types of buses* are P and O at load buses* e and f Cor V and 65 
at swing bus and P & V at generator buses, where at any bus 

*p» 


Vg » e^.f® 

P P P 


ca. 85 


At a P-<5 bus* since P and Q are sp>ecified, the mismatch 
vectors can be calculated as : 


pSp pcal 

p" P 


ca. Q5 




ca. 105 
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For P-V bu&> Bine® P and V are specified, bhe misma'tch 
vectors are calculated at any bus p using equation C2. dZ) and 
C3. IID given below 


AV? = 

P p p 


ca.iiD 


Where superscripts *sp* and *eal' stand for specified and 


calculated quantities resp>ecti vely. 


from equations Ca. 6Z> and C2. 71). 
ce. 85. 


.cal 


and are obtained 

P P 

is computed using equation 


Since load flow equations CS. 65 and C2. 75 are essentially 
non-linear algebraic equations, it can be solved by using any 
numerical technique. Solution of these equations using 

Newton-Paphson technique involves a linearized relationship 
between mismatch vector and voltage correction vector. For a 
n-bus system having bus-1 as slack bus, S to m as P-V buses and 
rest Cm+1 to n5 as load buses, the relationship between mismatch 
and correction vectors using NR method in rectangular 

co-ordinates, can be written as : 


« 
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C2. 125 


The shorty t.he ebcive equetion can be vn^meri as shown » 




where 
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ifi e£l.ll«rd 
correction 
fiubma.tr i ces 

CJg^pq 

CJ^^pq 

C-Tg^pq 

CJ^^pq 


T 

the Jacobian matrix and CAe» Atl Torm the 
vector. and are 

and their elements are defined as : 


as 

Far t^l fill 

derivatives 

of 

P V/. 
p 

r. t. 

e ' s. 

q 

= 

Par'Ll al 

derivatives 

of 

P w. 
p 

r . t. 

f 's. 

q 

= 

FarLlal 

der i vati ves 

of 

p 

w. r . 

t . « ^ s . 

q 

sr 

ParLlal 

deri vati vefi 

of 

P 

w. r . 

t. f /«. 

q 

s 

Partial 

derivatives 

of 

w. r . t. 

a •&. 

q 

sr 

Partial 

derivatives 

of 


r . t . 

f 'S. 

q 


vol tage 
Jacobi an 


It can be shown » that» the off-diagonal elements of Jacobian 
fiubmatrices Ci.e. q K p3 are as follows : 



as 

« G - f B 

p pq p pq 

C2. 143 

cJg^pq 

= 

« B + f G 

p pq p pq 

C2. 1B3 

C JgDpq 

= 

o 

C2. 163 

J^^pq 

as 

o 

C2. 173 

c^pq 

s 

f G + e B 

p pq p pq 

C2.183 

cjg^pq 

== 

f B — ft G 

p pq p pq 

C2. 103 


The diagonal elements of Jacobian fiubmatrices Ci.e. for q= p3 


are given afi» 



1 ' 


cjppp 


Se 6 
P PP 


n 

Z Ee G + f B 3 
q pq q pq 

q=l 

q/p 


C 2 . 20 D 


CJ :)pp = ^ 

•= P PP 


CJ 5 pp = a® 

c* p 


rt 

Z Cf G - e B 3 
q pq q pq 

q=l 

q/p 


C 2 . 21 J 


C 2 . 22 J 


C J^Jpp 

C jgjpp 


c JgJpp 


2 f 

P 


B - V cr G - e B 3 

p PP <=1 P^^i q pq 

q=l 

q^ 

n 

+ ) te G + f B 3 

P PP £__ q pq q pq 

q=l 

q? 1 p 


C2. 23J 

C2. 24J 


C2. 255 


The solution of equation C 2 . 135 gives us the Ae*s and Af’s 
which are used to up>date earlier estimates of e*s and f*s and the 
process is repeated till the mismatch AP and AO for P-0 buses and 

s 

AP & AV Tor P-V bufi«*s eonv«rrg<fe *to a pr^-assigned 'tolaranca valua 

j£?. 


S. 3, 1 Steps Tor Load Flow Solut^ion using Nfewton-Raphson Met^hod 


1 n Rec t angul ar Coor dl nates : 


The procedure for solving the load flow problem using Newton 
Paphson method in rectangular co-ordinates is described through 
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riow chart, given fig. C2. ID. Major compute ali onal sleps are as 

gi ven bel ow : 

1. Read system data including line data, loads at various buses, 
generator data. 

2. Form tY 3 matrix. 

BUS 

3. Initialize the voltage at each bus. If nothing is 
prespecified, assume a flat start Cl + JOD i.e. assign V at 
all buses c-qual to 1.0 p-u. and angle <5 equal to slack bus 
angle i.e. zero. Set iteration count k = O. 

4. Calculate P and Q using equations C2. 6D and C2. 7D with 

available values of e*B and f*s for P-0 buses, 
cal 

G. Calculate P , and voltage magnitude using eqns. 

C2. 6D and C2. 8D with avaialble values of e*s and f*s for P-V 
buses . 

5. Calculate px:iwer mismatch AP for all P-O and P-V buses and AO 
for all P-O buses using equations C2. Q2> and C2. lOD 
r espec t i vel y . 

cal 

7. For P-V buses, check if O ia within the sprecified O limits 
CO”^*' S i O”**^. If it is so. calculate A V ® using 

eqn. C2.11D. But if O**"^ is less than or greater than 

O » we proceed as follows : 


ir 

substitute 

if 


o'** 

> 

^max 



^max 


< 

^min 

a'** 

s 

q”*'’ 


» 


» 


subfili 'tu'ta 



START 


load flow da 


ta Sc f o nn Ty 


Assun^-^„ V .2...n 

^ slack bu c; ^ 

^ — "— ---- — - ^ 

'”■ gjt iteration onnr; -I- j^_.q 

— bus conni- _ 1 


Compute P^, 


P P 

apK pSP _ pk 
P P n 


PV 

type 


Compute 


bus p a PV or PQ 




( Vp ) 


pa 

^ t;^pe 
V "Computer” 
}!^ * QSP- 


Qmirr 


ComputeAQ^=:^ax_, qH 


compute 
( vSP )2 


Compute aQ^* 


(last bus! 


^0t0irrnin0 

inaxir>uin mi sn atch, 
max . [a pk^ ^ Q K 


Set p= p+l 


^ jgQ — -555il:^^smatch^ £ ? 

f Assemble Jacobian TjT” f — , 

’ X P 1 Compute and print line 

Determi ne F f l ows/ power loss,voltaqe 

compute Ae ^ and A f p ~ t, ^ 

I... i —■* ( STOP ) 

Compute - -X ^ . k 1 ^ 

^ ^ — 

Advance itn. count K=K 4-1 f 


«i??aSSSfa\^?oUg?SlSgj?iP^3on .ethod i 
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and t-ra-at, t.hls bua as a P— Q bus. Raact-lve power mismat.oh AQ is 
found using eqn. C2.105, 

8. If AP and AO at all P-Q buses and AP and AV® at all P-V 
buses are less than a pre-specif ied tolerence e, go to 
step-13, otherwise go to step-O. 

Q. Calculate the elements of Jacobian using equations CS. 

through C2. 255 and find the inverse of Jacobian matrix. 
lO Solve equation C2. 135 for voltage mismatch vectors Ae and Af. 

11. The assumed or calculated bus voltage components e and f are 
updated at all the buses using the correction vector 
obtained, as follows 


Ck+15 k . k 
e = e ‘4- Ae 


C2. 265 


^Ck+i5 ^k ^ .^k 

f = f + Af 


C2. 275 


k k 

vhere, e and f are the values of voltages at the kth 
iteration. 

12. Advance the iteration count by *1* Ck = k+15 and go to 

step-A. 

13. Load flow convergence obtained, calculate slack bus power, 
line flows and line losses etc. 


2. 4 Newton-Raphson Method in Polar Coordinates .- 


The load flow problem using Newton-Raphson method in p>olar 
coordinates CNRP5 can be formulated as follows. For any bus-p, we 



have 


2 


v“ = V 

p p p p 


V = V y 

q q pq 


lY le-J^P^ 
' pq* 


C2. 283 


where 6 Is Ihe angle of Ihe bus volt^age and 8pq is bus admi-tbanee 
angle. Then for any bus p, 

n 


■“p = ''p 


I 

q=l 


CY V 3 
pq q 


C2. 283 


Subs-ti t-utl ng relations C2. 283 in eqn. C2. 283 we get 
^ -jce + <5 ~ J 

] C2.303 


P P 


= I 


-JC8 +6 - <5 3 

C V |Y |V e ^ ^ 

pi pql q 


q=l 


Thus 


P = Real 
P 


n 




-jCd + S - 6 3 

^ pq p q 


q=l 


= I 

q=l 


Y |V Cos C& + S - <5 3] 
pi pql q pq p q 


n 


= I 


V V CG CosC<5 -6 3+B SinC«5 -6 33 C2. 313 

P q 


q=l 


q pq 

p q pq 

r 

1 r—' 

-jce 


and Op = Imaginary 


Z V |Y |V e 

pi pql q 


pq P q 


q=l 


= I 


C V V |Y I Since +6 - 6 33 

p ql pql pq p q 


q=l 
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n 

= ) V V [Q SinC6 -6 5-B CosC6 -6 5 3 

Z_^pqpq pqpq pq 

q=l 


C2. 325 


Equat/lon C2. 315 and C2. 325 form the pair of load flow 
equations at each bus in the system. The solution of load flow 
equations in pnolar coordinates involve real flower mismatch CAP5 at 
P-Q as well as at P-V buses and reactive power mismatch CA<25 at 
only P-Q buses. At any bus -p, these can be calculated as» 

AP = - p*^*^ 

P P 


AQ 


^sp ^eal 

V -s 


P and Q are obtained using eqns. C2. 315 and C2. 325. 
P P 

For calculated values of mismatch vector > the voltage corrections 
AV and A6 can be computed using NP relation in polar coordinates 
as shown . 


■ AP 


■•^1 

^2 ■ 

AQ 


/3 



Ad 

AV 


C2. 335 


where » 



J 

J 


2 

A- 


Is called the Jacobian matrix. » Jg* Jg and are Jacobian 

submatrices whose elements are <iefined as : 


C 5 p>q 
cJg^pq 


Partial derivatives of P^ w. r.t. 

P 


Partial derivatives of P^ w. r.t. 

P 


S 's. 

q 

V 's. 

q 
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CJ^^ppq - Pa.r'tla.l derl va.t.i ves of w. r.^. S^'s. 

CJ_,5pq - Partial derivatives of Q w. r.t. V 's. 

A p q 

For a n-bus system having Cm-15 P-V buses, the size of 
Jacobian matrix is C2n-m-l!> x C2n-m-lJ. 

It can be shown, that, the off-diagonal elements of Jacobian 
submatrices Ci.e. for q pj are. 



= V V tG SinC<5 -S J 

p q pq p q 

- B CosC6 -6 53 

pq p q 

C2. 345 

CJgJpq 

= V EG CosC6 -6 5 + 

p pq p q 

B SinC6 -6 5 3 

pq p q 

C2. 355 

CJgJpq 

= -V V EG CosC6 -6 5 

p q pq p q 

+ B SinC6 -6 5 3 

pq p q 

C2. 365 

CJ^Jpq 

= V EG SinCd -6 5- 

p pq p q 

B CosC6 -6 53 

pq p q 

C2. 375 


The diagonal elements of Jacobian submatrices Cl . e. for q = 
pD are 


CJ, Jpp = V V V CB CosCd -6 J-e SinCd -6 > 3 C2. 3SJ 

1*^*^ ^pqpq pqpq pq 

q=l 

q/p 

n 

CJ^3pp = 2V G +VvCS CosCd -6 3+B SinCd -d D 3 
2^^ ppp^qpq pqpq pq 

q=l 

q?^p . 


C2. 3d^ 
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V [G CcsC6 - <5 5+B SinCtS -S !>] CB. 402 
Pq p q pq P q 


CJ^:>pp 


2 V B + ) V 

p pp Z- q 

q=l 

q/p 


[G SinC6 -6 5-B Cost <5 -<5 5 ] 
pq p q pq p q 


C2. 41D 

Th* eolut^ion or oqn. C2. 33Z> givoa ua t,h» A^'s and 4V*b whii^ch 
are used t^o updale^ 'the earlier estimates or 6*s and Vs and iB'ie 
process is repeated till the mismatch AP for all P-V and P-Q bus«es 
and AO for all P-O buses become less than a pre-assiqned toleranece 
value 


2.4.1 Steps for Load Flow Solution using Newton-Raphson Mebhod din 
Pol ar Coor di nates 

Flow chart given in fig. C2. 2!> describes iha procedure test 
solving load flow problem using Newton-Raphson method in pyolar 
coordinates. However, major computational steps are as givcsn 
below : 

1. Read system data including line data, loads at various buses 
generator data. 

2. Form £Y 3 matrix. 

BUS 

3. Initialize the voltage ^ at each bus. If nothing Ms 

prespecified, assume a flat start i . e. assign V at all buses 
equal to 1.0 p. u. and angle equal to slack bus angle i .e. 


zero. Set iteration count k = O. 
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f'*'* *{ 

4. Calculate P and Q using equations C2- 31> and C£. 32Z> 

with available values of V*s and 6^s. 

B. For p-V buses replace temporarily by but leave <5^ as 

it is. Compute jTor p-V bus using 

cal 

©. For P-V busast check if Q is within the specified Q limits 

Ci.e. a: J£ If it is so, put of end 

1 

calculate AP only using eqn. C2.05. But if Q hits the O 
limit, proceed as follows : 


if 


> 

o”*"* 

put 


= 

Q»*X 

if 

o'*" 

< 

o""" 

put 

o'*" 

ss 

o”"" 


and treat this bus as P-Q bus. Use v as it is and find 
power mismatch AP and AQ using eqns. CS. and C2.10!> for all 
P-Q buses. 

7. If AP*s and AQ’s at all P-Q buses and AP*s for all P-V buses 
are less than a pre-specified tolerance go to step-12, 
otherwise go to step-8. 

8. Calculate the elements of Jacobian using eqns. C2. 34J through 
CS.AID and find the inverse of Jacobian matrix. 

9. Solve eqn. C2. 33J for voltage correction vector AV and Ad. 

10. The assumed or calculated voltage magnitudes and their angles 

i 

are updated using the correction vector obtained as follows : 


,Ck+lD 

O 


+ Ad 


yCk+lJ 


C2. A2::> 


V*' + AV^ 


C2. A3J 



27 


whAr*, S And V Ar« ■Lh* vaIuab aI. kt.h it-ArAt-ion. 

11. AdvAnce the Iteration count by *1* Ck = k+lD and go to 

step-4. 

12. Load flow convergence obtained, calculate slack bus p>ower » 
line flows and line losses. 

2. 5 Sample System Studies : 

The Newton -Raphson method outlined in section 2.3 and 2.4 
have been tested by solving the load flows for the following 

sample power systems using ECC-lOQO computer. 

1 . 1 EEE - 1 4 bus test system. 

2. IEEE - 30 bus test system. 

The specification of the above systems are given in 
Appendices A & B. 

A convergence critierion of 0.0001 p. u. for active and 
reactive piower mismatches and for voltage magnitude square 
mismatches has been taken on a 100 MVA base. The results have 
been obtained for both well behaved and ill-conditioned cases. 

1 1 1 -condi ti onl ng in the system has b€ren created by raising the 
series resistance of all the lines simultaneously from their base 
values dncreased R/OC ratio5. Load flow solutions on these 

systems using the flat voltage start are listed in Tables 2.1 to 
2.3. The number of iterations ‘and total C. P. U. time for load flow 
solution is shown in Table 2.1. Tables 2.1 and 2.2 pier tain to 
unadjusted solution Cno O-limits imjsosed^. whereas Table 2.3 
corresponds to adjusted solutions considering the generator 
0-limits. Table 2.2 compares the convergence characteristics for 
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i 1 1 “Cioridl li oned Ei-LuaLion. Cases for increased R/X raLio of lines 
have been shown as CC1 +<dDR + jX3 in Table 2. S which indicates that 
the series resistances of all the lines are increased by times 

thei r base val ues . 

The convergence pattern of both NRR and NRP methods have been 
studied by plotting the values of power mismatches obtained in 
each iteration for the study of unadjusted base case and are shown 
in figures C2.3J and C2. Analysis of these figures clearly 

reveal the quadratic convergence pattern of load flow solution 
using either NRR or NRP method. It is also seen that with the 
increase in number of buses* the NRP seems to be better than that 
of NRR as for as CPU time is concerned. However the number of 
iteration required in both these methods are same and do not vary 
much with the increase in problem size and kind. 


TABLE - 2.1 


Comparison of NRR and NRP methods for load flow solution 

C unadj usted case!) 



IEEE - 

14 BUS SYSTEM 

IEEE - 

30 BUS SYSTEM 

Hfirt^hod 

No. or 

I 'tor at ion 

Total C. P.U. 

Time C seconds!) 

No. or 

I tor at.ion 

Total C. P. U. 

Time (! seconds!) 

NRR 

3 

2.085 ‘ 

3 

15.564 

NRP 

3 

2. 631 

3 

18.071 



1.00 



FIG - 2.3 Convergence of Newton Rophson 
method in rectangular Sc polar co-ordinates 
(IEEE — 14 Bus System) 




oon 



FIG - 2.4 Convergence of Newton Rophson 
method in rectangular Polar co-ordinates 
(IEEE - 30 Bus System) 
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2. e Discussion and Conclusion : 


From the result shown in above tables » it is clear that NRR 
is slightly faster than NRP. Number of iterations required in 
both cases are almost the same. However number of equations and 
variables in NR method in rectangular coordinates is greater than 
that in NRP, by the number of P-V buses. Hence memory requirement 
in NRR is more than that in NRP. Comi.?utati on time in rectangular 
coordinates CNRRD is less because time consuming trigonometrical 
functions are not present. Thus, the formulation in rectangular 
coordinates takes less computational effort and that in 
polar -coordinates takes less memory space. Though the number of 
iterations required to obtain load flow solution using either of 
the methods is almost independent of problem size and kind. 

Newton *s methods are also unaffected by ill-conditioned 
situations viz. negative transfer reactances which can not be 
netted with positive reactance of the line and a situation where 
high and low impedance branch terminate at the same node. The NR 
method needs 4 to 5 iterations to reach an acceptable solution for 
a large system. However, the time per iteration in both these 
methods CNRR and NRPI) increase almost directly as the number of 
buses of the network. If Q-limits are considered, the NR method 
takes more iterations than that for unadjusted solution. 

Since, the NR method requires number of iterations 
independent of problem size and kind, it is more advantageous for 
large scale systems. NR method has quadratic convergence 
characteristics and is the best among all methods from the stand 
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point of convergence criterion. With NR method, the elements of 
the Jacobian are to be computed in each iteration, and hence, the 
time per iteration is considerably longer than Gauss-Seidel method 
C about 7 timesi). In fact, the NR method is more reliable and is 
able to solve a wider variety of problems. Its main drawback of 
the large requirement of computer memory can be over-come by a 
compact storage scheme. It can be made faster by adopting the 
scheme of optimally ordering of buses etc. Further in NR method 
it is necessary that the initial guess is not far off from the 
nominal solution otherwise there will be quite a large truncation 
error . 


i 




CHAPTER - III 


FAST DECOUPLED LOAD FLOW METHODS 


3. 1 Introduction : 

Amongst the various load Tlow methods avail able» the fast 
decoupled load flow method is the most popular with utilities for 
obtaining load flow solution because of its reduced memory 
requirements and providing faster solution. The FDLF method is 
attractive for accurate as well as for approximate off- and 
on-line routine and contingency calculations for networks of any 
size and kind. In this chapter » the fast decoupled load flow 
CFDLFZJ method of Stott and Alsoc CIS] in polar coordinates and 
existing FDLF method in rectangular coordinates [21,32] are 
discussed. New FDLF algorithms in rectangular coordinates 

considering certain modifications in the existing one is also 
presented in this chapter, which is found to provide convergence 
especially for ill-conditioned situations in the system. Various 
existing and new models of FDLF method have been tested on 
I . E. E. E. 14-bus and 30-bus systems and their performances compared 
for well-behaved and ill-conditioned cases in the systems. 

3. 2 Newton-Raphson Decoupled Load Flow Model 

Newtrori— Raphson load flow met<hod developed by Tlnriey and Hart, 
C63at B. P. A. Portland CORE-USA3 is one of the most powerful 
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m 9 t.ho<i 8 . Ther number of i'tora.'tions ro<]uir 9 ci ■bo ob'La.in soIu'Lion by 
bhis method is independent of problem size and kind. Even with 
flat voltage start, it takes nearly 4 to 5 iterations to obtain 
the solution, even though the time taken to coital ete one iteration 
is nearly seven times to that of the Gauss-Seidel method, still 
this method is faster for large scale systems. This method is not 
affected by ill-conditioned situations and also with slight 
modifications, all the other factors such as representation of 
transformers etc. can be includcrd in the same Jacobian. 

The load flow equations used in this method for computation 
of voltage corrections are given as. 


AP 

AQ 


where. 


H 


km 


J 


km 


H 


N 


AA 

AV 

V 




N, 


km 


V 

k m 
^m 


m 


km 




m 


9W 


m 


C3.1D 


H, J and L are the submatrices of the Jacobian. 

The decoupled method Cl 43 is actually an extension of 
Newton method in j>olar coordinates €Kq[n. C3. ID with certain 
appr oxi mat i ons . Si nee the changes in r eal power C i , e . APD is 1 ess 
sensitive to the change in voltage magnitude Ci.e. AVD and 
changes in reactive power AQ is less sensitive to the changes in 
angle A<5ii we can write the eqn. C3. ID in the following form. 
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AP 


H O 


A<5 





AV 

AQ 


O L 



. 




V 


C3. 25 


EquaLioris C3. 25 are decoupled equate! ons which can be expanded 
as 


CAP3 


CAQ] 


CH3 


tLJ 


[ A<51 

AV 

V 


C3. 35 


c:3. 45 


li can be shown t,hat, the orr— diagonal elemen'ts of Jacobian 
C i . e . f or q p5 ar e • 


H = L = V V CG SinCA -<5 5-B CosCA -6 5 3 
pq pq p q pq p q pq p q 


C3. S5 


The diagonal eletnent/S of [H3 and [L3 are 



B 

P PP 



C3. 65 


PP 




B 

P PP 


+ Q 


C3. 75 


Equations C3-35 and C3. 45 are formulated and solved 
successively. The latest values of S found from eqn. C3.35 are 
used to solve for eqn. C3. 45. The decoupled method converges as 
reliably as the full Newton -Raphson method, although it takes more 
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number of iterations as compared to the Newton's method. The 
saving in terms of memory requirements is about 7S>i for Jacobian 
elements although overall saving of the memory is only of the 
order of 40-50>S. The computation time per iteration is also 
10-205< less than full Newton -Raphson method. 

3.3 Fast Decoupled Load FI ow Method in Polar Co-ordinates : 


The decoupled method [14] can be further simplified without 


appreciable loss of accuracy. In practical p>ower system, the 
following assun^tions hold good : 


1. C6 -6 D is small. 

P q 

2. G SinCd - <5 J << B . 

pq P q pq 

3. Q << B .V®. 

^ PP P 


Applying these assunqations to equations C3.3^ and C3. 4I>» we 


get 


CAP] = CV. B'.V] CAA] C3.SJ 


CAQ] 


CV. B". V ] 


AV 

“V“ 


C3. 03 


The elements of the matrices CB'] and CB''] used in eqns. 
C3. 83 and C3. 03 are directly derived from elements of E -B] . The 
decoupling process is given a 'final form by following additional 
assumptions : 

a3 Omitting from CB'] the representation of those network 
elements that affects MVAR flows i . e. shunt reactances and 
off-nominal in-phase transformer taps. 
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b5 Omi'tting from the angle shif'ting effects of phase 

shifters. 

c? Neglecting series resistances in calculating the elements of 
tB'] and dividing eqns. C3.83 and C3. 95 by vdien 

calculating for pth bus, and setting = 1 p. u. 

With these assumptions, equations C3. 85 and C3. 95 take the 
following form 

CB*3 1662 C3.105 

EB**3 CAV3 C3.115 

vdiich are the set of equations used in fast decoupled load flow 
model C15] in polar co-ordinates. 

In eqns. C3. 105 and C3.115 CB'] and CB"3 matrices are real, 
sparse and have the structure of [H3 and CL3 respectively. Since 
they contain only network admittances, they are constant and need 
to be triangulated once only at the beginning of the iteration. 
Using the constant triangular factors of tB’3 and CB”3, very fast 
repeat solutions for [Ad3 and [AV3 are obtained. The best scheme 
for F1X.P is to solve eqns. C3.105 and C3. 115 alternately, always 
using the most recent voltage values. Each iteration cycle 
coR^rises one solution for EA53 to update CA3 and then one 
solution for [AV3 to update [V3, and termed as Cl~6, 1-V5 scheme. 
The iterative procedure is repeated till the real and reactive 
power mismatches become less than a pre-assigned tolerance value 


■ AP 
"V“ 


AQ 

“V 
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Varloufi sLorage schcrmas can be applied to take the advantage 
of sparsity and symmetry of CB*3 and CB*'3. However* CB'] is 
symmetrical only if the phase-shifters are neglected. 

3, 3. 1 Steps for the FI3LP iteration scheme : 

Major computational steps for solving load flow problem using 
FIX-P method are shown in a flow chart Cfig. 3.15 which are 
summarized as follows : 

15 Read system data including line data* loads at various buses* 
generation data etc. 

25 Form Y-bus matrix* CB*3 and tB'*] matrices. 

35 Initialize the bus voltage to Cl+j05 p. u. for all P-0 buses 
cj) 

and CV + J05 for all P— V type buses. 

45 Set Iteration count k = O. 

55 Calculate bus real power P and mismatch AP for all buses 
except slack bus using 

AP = P®^ - P’^*^ C3.125 

65 Check* if real piower has converged to its sF»eelfied values 
CMax. AP < /f5. If so* go to step CS5* otherwise go to step-7. 
75 Advance iteration count by O. S Ck = k +0.55 and solve for A<5 
using eqn. C3.105. 

S5 Update the voltage angle ,«5 at all buses C except slack bus5 
using the corrections A«5 obtained in step-7 

& = 6 + Ad C3. 135 

S5 Calculate bus reactive power O and mismatch AO for all P-O 
buses using following eqn. 



Compute inverse if C®.!! ^ L® J 
and ifeit-ialize bus voltage 


Initialiie Itn.. count K = 0,0 


Compute vectors P and ^P 
for all Duses except slack 



Compute vector v 

AV from Eqn. [Advance K = K -f o«5 

(3.11) and ,j 

update V « V Compute vectors AS from £qn 3.10 

“ I and update S * 6 + S 


Compute vector Q and 
AQ for P-Q type buses only. 



Advance 

K « K + 0.5 


A Q max I ^ ^ 7 


Convergence obtained 
Compute and print 
desired L.F. results 


[a P max 


STOP 


Fig- 3.1 Flow chart for FDLP Algorithm for Unadjusted 
case using 16-lV Scheme. 













105 Check, if reactive powers at P-Q buses have converged CMax. 

AO < c5 . If so, go to step-14, otherwise go to step-ll. 

115 Advance iteration count by 0.5 Ck = k + 0.55 and confute 
vector [AV] from eqn. C3. 115. 

125 UfKdate the voltage for all P-Q buses using the corrections 
AV obtained in step-11 as 

V = V + AV C3. 155 

135 Go to step-5 and repeat the procedure. 

145 Check, if max. AP < £. If not, then go to step-5 and repeat 
the proecedure, otherwise go to step-15. 

155 Load flow convergence obtained, calculate slack bus power, 
line flows & line losses. 

3. 4 Fast Ctecoupled Load Flow Method in Rectangular Coordinates : 

Fast decoupled load flow method Cl 53 in polar coordinates 
CFEX-P5 discussed earlier is found to be the most popular method 
with utilities because of their lesser memory and CPU time 
requirement. The formulation of a load flow method in rectangular 
CO— ordinaters, as observed in Chapter— 2, offers faster solution as 
compared to its formulation in polar co— ordinatefS , mainly because 
of the reasons that it does ncit involve additional computational 
efforts in solving trigonometric functions. With this in view, 
development^ of fast decoupled load flow methods in rectangular 
co-ordinat,l|rs CFiaLR5 appears to be more attractive. Literature 
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survey shows ‘Lh&'L s ll'L'Lle effort, h&s heen ina.eie in developing 
FDLR models. The two efforts one by Arumugam [213 and other by 
Singh et al . C22] ineidently result into similar formulation of 

FDl-R model . The development of FDtR model as suggested by them» 
using two different approaches, one from FDLP model using linear 
transformation and the other from Newton-Raphson load flow 
equations are described below. 

3. A. 1 Derivation of FraJR from F13LP by linear transformation : 


The real and Imaginary parts of the pth bus voltage in terms 
of its magnitude and angle are 


e = V Cos <5 ; f = V Sin 6 C3. 165 

p p p • p p p 


By taking the exact differentiation of the equation C3. 165, 
we get 


Ae 


Af 


-V Sin S Cos 6 


AS 

P P P 


P 

V Cos S Sin A 


AV 

P P P . 


P 


C3. 175 


Inverting eqn. C3. 175, we obtain 


- 


n S 

Cos 6 



P 

P 

AS 


V 

V 

P 


P 

P 

AV 


Cos S 

Sin 6 

P 


P 

P 


Ae 

P 



C3. 1S5 
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Applying the transformation C3. 1S> in eqns. C3. 105 and C3. 11> 
in rectangular coordinates the following equations are obtained : 


■ AP 

“V" 

= CB'] CA^ 


Afir 

m : 

C3. 105 

t 1 

= [B"] [A!' 

1 

^23 

Ae 

Af 

C3.205 


where CA’3, ^^2^’ and tAg] are diagonal matrices with 

ppth element being given as 


ca;5 

1 pp 


Sin «5 

^P 

V 

P 


CA1> 

a pp 


Cos <5 


ca:'> 

1 pp 


Cos <5 


and 


2 pp 


Sin 6 


In the polar version, the voltage magnitude equation do not 
appear for the P“V buseSi Hence, the equations C3. 195 and C3. 205 
in rectangular coordinates have to be augmented by 

[A V®] = tA^'» A^’3 


■ Ae ■ 

W 

L 


C3.215 
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Where, CA”*3 and ^^“*3 are diagonal ma'trices wilh ppl-h 
elemenl given by 


CA"*:> 

1 pp 


On combining equations (13. 102), <13. 30^ & C3. SID and making 

the same assumptions as made in FEA-P, namely V ssl.O, 6 -6 m O 

P P <5 

etc. and after some manipulations, vas get 
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22 2 2m m 2,m+l m+l*^ 2n n 
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nm m 


C-B^ C-B* r 3 

ri*m4*l m-**l nn n 




CB* * :> 
2m m 


2 > m+1 m+1 


CB' * 3 
2ri n 


^®A2V 


CB' * 5 
nm m 


CB' « :> . 

n»m+l n+1 


CB' « > 
nn n 


m+1 , m+1 m+1 !) 


CB*' , * 5 

m+1 » n n 


0 CB" e ^.3 CB" e D 

n,m+l m+1 nn n 




o CB" ^ r .3 CB" r 5 

n » m+1 m+1 nn n 


r sf ^ 

‘0 

O 

o 


• 

• 

- 

o 

af 

0 

b 
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In t/h® Ab>ov® d®ri VAt-i on » 'Lhi® bus~l Is t>&k®n &s slAck busr 
buses a-m on P-V buses end buses m+1 to n es P-Q buses. To 
echieve decoupllngt the following further essumptlons ere nede : 

ID A flet voltege stert Ce+Jf = l.O+JO. OD is used. 
ilD The met rices obteined using the flet voltege stert ere 
treeted es constant et every subsequent iterative step. This 
avoids the re-ealeuletion and re-tri angulation of the 
Jacobian metric at each iteration. The algorithm obtained, 
therefore, utilizes the fixed tangent slopes corresponding to 
the flat voltage start instead of the varying tangent slopes 
at each iteration. 

With the above assumptions, eqn. C3. aSD is reduced to the 
following set of decoupled equations in rectangular coordinates : 



= CB*3 

^a 

ca.asD 

AP 

n 


Af 

n 



= [B'‘3 


C3. BAD 



A«F 

n 


1 


Afip 

2 ! 


AV^ 

m 

= CB"'I 

A«r 

m 

C3. aSD 


= atI3. 


where t B*' * 3 



•» / 


Since tAf] and hence Cfl are known for the P-V buses from 
eqn. C3.e33. £e3 for the P-V buses can be calculated using its 

specified voltage magnitude from 


[f^i 

Replacing eqn. C3. 2S:> by eqn. C3. 2e:>, 
equations become. 


C 3.265 

the final FISLR 


CAP] = [B‘] CAfl 


C3. 275 


C Ae] 


C C V^^] ^ 




C3. 285 


CAO] = CB"] [Ael 


C3. 2S5 


The above equations can also be obtained from Newton-Raphson 
load flow equations in rectangular coordinates as indicated below. 


3. 4. 2 Derivation of FCl^R from Newton-Raphson load flow equations : 


With the flat voltage Cl.O+jO. 05 start at all the buses and 
also using it for forming the Jacobian in NR load flow equations 
C2. 135, its elements get modified as follows : 


OP 

P 


q 


% 


dP 

P 


Of 


q 


% 

q 


C3. 305 


C3. 315 
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OP 

P 

P 


& + G 

P PP 


C3. 325 


OP 

P 

P 




% 

■S? — ~ A — G 

Sf p PP 


where > 


and 



C3. 335 


C3. 34,5 


C3. 355 


C3. 3S5 


C3. 375 


From equat.lonB £3. 365 and C3. 375, we may assume a << 6 and 

P PP 

b << B and hence a and b can be neglected in equations 

p pq p p 

C3. 32-3. 355. Further, neglecting G and 6 in the above 

pq PP 

expressions, the load flow equations reduce to. 


CAP3 = CBp,] [Af] 

tAQ] = CB^] [Ae] 

C AV^3 = 


C3. 385 

C3. 3Q5 


CB*'n [Ae3 


C3. 4,05 
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Where, tBp] and t,he elements of C -B3 matrix. 

Incorporating various assumptions of FDLP in forming matrices CBpl 
and ‘V- they get modified to matrices [B*3 and tB'*3 respectively 
cpf FESLP model of Stott and Alsac C1S3. Replacing CBp3 and CB^3 by 
[B'3 and CB"! matrices respectively in eqns. C3.38D and C3. SBD 
and using eqn. C3. 265 in place of eqn. C3. 405, we get the final 
load flow equations for FDLP model which are nothing but those 
obtained from FCA-P by linear transformation reprcfsented by 
equations C3. 27—3. 285 . 

3. 4. 3 Iterative Scheme for Fast Decoupled Load Flow in 
Rec t angul ar Coor di nates 

Equations C3. 27-3. 285 are solved successively using the most 
recent values. The best Iterative scheme found to be Clf, le5 
scheme, which comprises of one solution for [Af3 using equation 
C3.275 to up>date Cf] for both P-V and P-Q buses followed by one 
solution for [el using equation C3. 2S5 for only the P-V buses and 
finally one solution for tAel using eqn. C3.;^5 to up>date Ce3 for 
only the P-Q buses. Iterations are repeated till the real and 
reactive p>ower mismatches are within a specified convergence 
tolerance. Major computational steps for solution of load flow 
problem using FDLR model are given in a flow chart Cfig. 3.25 
which are summarized as follows : - 

15 Read system data including line data, loads at various buses, 

gen^feratlori data ©tc. 






Fig- 3.2 Flow charffor FDLR Algorithm for unadjusted 
case using ( If - 1®) Scheme. 
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25 

35 

45 

65 

65 

75 

S5 

Q5 

105 

115 

125 

135 

145 

155 

165 


Form y-buB matrix, tB'J and IB”3 matrleas. 

Initialize the bus voltage to C1.0+J0.05 p. u.for all P-Q and 
P-V buses. 

Set iteration count k = O. 

Calculate bus real power P and mismatch AP for all buses 
except slack bus using eqn. C3. 125 

Check if bus real powers have converged CMax. AP < ^ 5 . If so, 
go to step-11, otherwise go to step-7. 

Advance iteration count by 0. 6 Ck = k + 0.65 and solve for Af 
using eqn. C3. 275. 

Update f at all buses Cboth P-V and P-Q buses5 using the eqn. 

f = f + Af C3. 415 

Calculte e f or all P-V buses using eqn. C3. 2S5 
Calculate bus reactive power Q and mismatch AO for all P-Q 
buses using eqn. C3. 145. 

Check, if reactive px>wers have converged CMax. AQ ^ *5. If 
so^go to Btep-16, otherwise go to Btep-12. 

Advance iteration count by 0.5 Ck = k + 0.55 and compute 
vector Ae using eqn. C3. 205. 

Update e at all P-0 buses using eqn. 

e = e -f- Ae C3. 425 

<3o to step-5 and repeat the procedure. 

Check, if max. AP < «. If not, go to step-S and repeat 
the proecedure, otherwise go to step-16. 

Load flow convergence obtained, calculate slack bus p>ower , 
line flows, line losses etc. .. lK' 

r;'";.:' ioioU.. 
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3 ® N»w Formulations of Fas! Decoupled Load Flow Method in 


Rectangular Co-ordi nates : 


Since, the e>d. sting FDLR model [31,323 described in earlier 
section considers all the assumptions of FISLP and in addition, 
certain further unrealistic simplifications, it is found to offer 
a convergence pattern inferior to the FDLP model. In fact, the 
number of iterations taken by FDt-R is found to be always more than 
FCX-P. Time superiority of FCA,R may be due to less time required 
in power calculations in rectangular coordinates than in polar 
coordinates. To overcome the drawback of the above model, two new 
formulations of FDLR model are developed as below : 

3.5.1 New FESLR Model -I CNFDR-I3 : 


By considering no further assumptions in the model to those 
already considered in FEA-P, for the pth bus, voltage in 
rectangular coordinates can be written as, 

^ s V C«o& & C 3 . 43^ 

P P P 

f = V Sin 6 
P P P 

Using exact differentiation, Ae^ and Af^ can be written as 

Ae = - V Sin 6 ^ A<5 + Cos 6 AV 

p p P P P P 


C3. A53 



And 


C3. 46:5 


Af 


V Cos S A<5_ + Sin S AV 
P P P P P 


Solving these eqns. , we obtain 


AV = Cos S Ae + Si n A Af 
P P P P P 


C3. 475 


and 


A<S 


Si n S Cos A 

P P 

V “V 


C3. 435 


Substituting values of Sin A and Cos A from eqs. C3. 435 and 

P P 


C3. 445 we get 


AV 


e f 

P P 

ry— Ae + ry— Af 

P P 


C3. 4Q5 


AA 


f 

P 

A* 

V^ 

P 


e 

P 

_ Af 
V^ P 
P 


C3.S05 


Writing FEX-P eqn. for real power mismatch for bus-p» we have 
AP 


n 


I 

q=2 


AA 

pq q 


<: * 1 » bei ng the si ack bus5 


or 


AP 


n 


q=2 


n 


f 

q 


n 

e 

q 


B* 


+ ) B* 


L 

q=2 

pq 

V® 

q 

q 

Z. M 

q=2 

V^ 

q 


r 

q 

-] 


n 

e 

q 

B‘ 

pq 

q 

q 


) B* C 

q=2 

V® 

q 


Af 


Af 5 

q 


C3. 615 
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Eqn. C3. BID nua.y b« wrlt,t»n in matrix form as 


CAP™] 

P 


CB* 3 CAf'“3 

pq q 


C3. B2D 


Similarly, the FE)LP equation for reactive power mismatch at bus-p 


1 . e. 


= I 


B" AV 

pq q 


Cfor P-Q buses onlyD 


q=m+l 


can be modified by substituting the values of AV^ from eqn. t 3 . 4 dD 


I “p 


q=m+l 


q q 

1" Ae + ^ Af D 

pq Vq q Vq q 


^ B" ^ Af I = V B" TJ- ^ 

qsm+l ^ J q=m+l ^ 


C3.53D 


Eqn. C 3 .S 3 D can be written in matrix form as 


_ - - Wl- 

CAO 3 
P 


CB" 3 tAe 3 

pq q 


C 3 . BAD 


Hence, final set of equations for NFDR-I including voltage 
irI SYRa'Lch i&quB.'Li onss for p— V bus(&s can wr i 't'L^ri as 



C3.55D 


CAP*”] = CB'l CAf*"] 

C.] = Ctv“P]® - 


CAQ*"] = CB"3 CAe*”] 


Where, the superscript m stands for modified values 


AP 


m 


AP 

P 

77- + f . C Ae:> 

V 1 

P 


AO 


m 


% 


— + CAf 3 

P 


Af 


m 


e 

q 

q 


A<er 


m 


=7— A& 

V q 

q 


Here» f^CAei and f^CAfS denote functions of 
1 3 

respectively* given as 


n 


f^ CAe3 
1 


= I 


, — 5 — Ae 

pq V® q 

q=2 q 


n 


fgCifJ 


= -E 


B” Af 

_ pq V ^ 

q=m+l q 


C3.S6D 


C3. S75 


given by 


C3.SS3 


C3. sci:> 


C3. €505 


C3.ei5 


Ae * s and Af * s 


and 



56 


During the first half iteration, Aa is not available, 
therefore, AP is taken as AP/V, and in eubeequent iterations when 
CAe3 and CAfl assume certain values, these are used to compute 
CAP 3 and C AQ 3. Kqna. C3. BB-3. B75 are solved successively using 
the most recent values. The best iteration strategy is termed as 
the Clf, le:> scheme, which comprises of one solution for tAf3 to 
update Cf] for both P-V & p-Q buses, one solution for [e3 using 
eqn. C3.^5 for only P-V buses and then one solution for CAel to 
up>date [e3 for only the P— O buses. Iterations are repeated until 
the real and reactive power mismatches are within a sp^ecified 
convergence tolerance. Major computational steps for the proposed 
model of FM..R is presented in a flow chart CFig. 3. 3:> which is 
summarized as following : 

15 Read system data including line data, loads at P-Q buses, 
generator data etc. 

25 Form Y-bus matrix, [B*3 and tB”3 matrices. 

35 Initialize the bus voltage to C1+J05 p. u. for all P“Q and 
P-V buses. 

45 Set Iteration count k = O. 

B5 Calculate C AP"*3 for all buses except slack bus using eqn. 
C3. 125. 

65 Test for convergence of bus real powers Cnax. AP ^ ^5. If 
convergence obtained, go to step-ll, otherwise go to step— 7. 
75 Advance Iteration count by O. 5 Ck = k + O. S5 and solve for Af 
using eqn. C3. 605. 

85 Upjdate f at all buses Cboth P-Q and P-V buses5 using the eqn. 

f = f + Af 



START ) 


Compute inverse of b' and b" 
Initialize bus voltage 




Ini tialize itn. count K = 0.0 


Compute vector P and A P for 
all buses except slack 


Compute vector 
from Eqs. (3,57) 
and (3,61) and 
update e * e +Ae 


Comnu te 


— ~^Mo Substi- 
K > tute 

^ Af = ( 


Advance 
K =K + 0.5 


K ^0.5 t^ubs ti tu te *0 


Compute 


Compute vectors Af from 
Eqs, (3.55) and (3.60) k 
update f * f +Af 


Compute vector e for 
PV buses using eqn.(3.56) 


Compute vector Q and 
AQ for P-Q type buses only 




Convergence obtainea 
Compute and print 
L.F, resu l ts 

r STOP ') 


A PmaxlA^^ 


Fig- 3.3 Flow chart for NFDR Algorithms for unadjusted 
case using If - le Scheme. 
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05 CAlculat* « for *11 p-V bucoa ueing aqn. C3.B65 

105 Calculat® AQ for *11 P-O buafta ueing *qn. C3. 575. 

115 Test for convergence of bus reactive powers CMax. AQ < 4 ? 5 . If 
convergence obt*lned. go to atep-lS, otherwise go to atep-12. 

125 Adv*nce iteration count by 0.5 Ck * k +0.55 and compute 
Ae using eqn. C3. 615. 

135 Update e at all P-Q buses using eqn. 

e = e + Ae 

145 Go to step-B. 

155 Check, if max. AP £ jc. If not go to step-5, otherwise to go 
step -16. 

165 Convergence obtained. Calculate slack bus power, line flows, 
line losses etc. 

3.5.2 Modified version of New FOLK Model -I CNFDR-IM5 


In NFE5R-I model described in section 3.5.1, and CB"] 

matrices are formed using all the assumptions made by Stott and 
Alsac in their FDLP C15] model. Van Amerongen C243 in his work 
has suggested certain modifications in usual FE4-P model which 
results into better convergence pattern esp>eclally for solving the 
ill-conditioned systems. The modifications suggested by him is 
mainly on considering a difi^erent combination of additional 
assumptions other than made by Stott & Alsac £153 in forming EB 3 
and CB"3 matrices. Hence a modified version of New FOLK model was 



5 


also t-ri^d out, ineorF»orat.irig t,h* changm auggABtod by Van 
Amarongan [343 whlla forming [B‘3 and [B*'3 matrieaa. Various 
modifications includa : 

a3 Line series resistance should be neglected in forming CB"3 
matrix instead of [B*3 matrix used in usual FDLP model C153. 
bD tB"3 matrix is formed considering double the nominal value of 
line shunt capacitances. 

The iterative scheme for solving io.»d flow problem using the 
modified version CNFDR-IM:> la exactly same as that for NFDR-I 
scheme described in section 3.5.1. 


3.5.3 New FDLR Model -I I C NFDR-I 1 3 


A different version of FDLR was developed which is directly 
derived from Newton -Raphson load flow equations in rectangular 
coordinates using the only approximation that, the initial value 
of voltage at all buses are assumed to be for calculation of 

elements of Jacobian. The elements of Jacobian will be as 
foil ows : 


1 pp 


1 pq 


3 pp 


Z 6 + G =: G 

pq pp pp 

q=l 


G 

pq 



B 

pp 


B . 
shp 


C3.S25 


C3. 63J 


C3.64) 


6 


CJ^:> = B 

2 pq pq 


» 2* 

2 PP P 


* O 

3 pq 


cjp = ar 

-4 pp p 


^ pq 


n 


cjp 

B pp 


B + > B 
PP ^ pq 
q=l 


B + B . 

pp Bhp 


S pq 


= B 


pq 


n 


cj^:> 

o pp 


= I 


(3 — G 

pq pp 


= - G 


pp 


q=l 


e pq 


--e 


pq 


Using abov«» the P-f eqn. Cfor P-0 and P-V busesD will be 


AP 


A A 

= I 


G Ae -I- CB 

pq q pp 


B . >Af + 

shp p 


q=a 


n 

I 

q=a 

h 


B AT 

pq q 


n 


or CAP 


- V' O Ae3=CB -B.^Af + ; B At 
p y F»q q pp el'P P ^ Pq q 

q=2 q=2 

4p 


C3. eB5 


C3. 665 


C3. 675 


C3. 685 


C3. 605 


C3. 705 


C3. 715 


C3. 725 


C3. 735 



or 


[AP"1 


CB'l CAT] 

r 


C3. 745 


where, the diagonal elements of [B*] = B 


PP ®shp ^ 


off-diagonal elements of [B* ] = b 

r pq' 


For the voltage update equation for P-V buses, we 


use 


Cel = CCV®^]® - 


C3. 765 


Similarly, Q-e eqn. for P-Q buses can be derived as. 


n 


AO = CB + B Ae-Vc Af 

p PP shp p pq q pq q 

q=2 q=2 

5^P 


or C AO 


Z 


n 


<3 Af ]= CB + B . 5Ae + ) B Ae 

pq q PP shp p / ■ F»q q 

q=3 q=S 


or 


c aq"*] 


[B"D CAe] 
r 


C3. 765 


where, the diagonal elements of = ^pp'^^shp off-diagonal 

elements of B" = B . 

r pq 

B Is the total additional shunts at bus p. Thus the final 
shp 

set of eqns. for NFE*-II model will be 


CB'3 tAf] 
r 


CAP*®! 


C3.775 
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EAQ®] * [B"3 CA*3 

r 


AP“ = 

AP + 

f uncti on 

of 

CAe5 

AQ 

AO + 

f uneti on 

of 

CAf5 


C3. 7S> 


C3. 705 


During t.h«r Tirsst h«.lf Itftr&tion, [Ae3 is not svailable. 
'thAr*ror«. AP™ is t,akttn aa AP, but, in eubaaquant, itarationa whan 
CAa3 and [Ai"3 assume certain values » these are used to coirqpute 
CAP'*^3 and tA0"^3. For load Tlow aolutlonClf, la5 schema described 
for NFC*:-I model in section 3.5.1 is also used for the pr opposed 
NFC^-II model . Iterations are repieated until the real and reactive 
piower mismatches are within a specified convergence tolerance. 

3. 6 Sample ^stem Studies 


The FDLP, FEH-R, NFE«-I , NFDR-IM and NFDR-II models discussed 
in previous sections have been tested by solving the load flows 
for the following jsower system examples on DEC-1 000 system. 

1 . I ESE 1 A bus test system. , 

2. IEEE 30 bus test system. 

The specifications of the above systems are given in 


Appsendices - A & B. 



Mat/ Hod 


IEEE 


14 BUS 


IEEE 


30 BUS 



No . o jf 

I tafra'tions 

CPU Time 

C seconds!) 

No. of 
Iterations 

Fca-p 

4. O 

0. 302 

3. 5 

FiaLR 

5. 5 

O. 410 

6.5 

NFE«-I 

11 . 6 

0. 615 

14.5 

NFE«-IM 

12. 6 

O. 628 

17. 5 


NFC«-II 


22. 5 


O. 856 


33. O 











TABLE - 3.2 


Efffect of Large Line R/^X Ratio; on Convergence CTypfe-a Simulations 


AS IEEE - 14 Bus System : 


111 -condi ti oni r»g 

oaae 

Cl+aSR + JX 

Humbi^r of I i«pr atrl ona 

FDLP 

FEH-R 

NFDR-I 

NFDR-I M 

NFI»-II 

l.B R + jX 

B. B 

8. B 

13. 0 

13. S 

Div. 

2.0 R + JX 

S. 0 

13.0 

14. 0 

14. S 

Div. 

2. 5 R + JX 

11 . O 

28. B 

16. 0 

16.5 

Div. 

3.0 R + JX 

14. B 

30.0 

10. 0 

20.5 

Div. 

3.B R + JX 

21 . O 

Osc. 

27.0 

28.5 

Div. 

4.0 R + JX 

3G. B 

Div. 

13B.0 

Div. 

Div. 


Ose . - Osei 1 1 ati ng 
Div. - Diverged 


BS I EEE — 30 Bus System : 


Ill -condi ti oni ng 

case 

C1+«SR + JX 

Number of I'torationa 

FEX-P 

FDLR 

NFDR-I 

NFE«-IM 

NFDR-I I 

l.B R + JX 

6. O 

11.5 

16. 5 

10.5 

Div. 

2.0 R -I- JX 

8. B 

IB. O 

18.0 

24.5 

Div. 

2.B R + JX 

12. 0 

42.0 

20. 0 

38. 5 

Div. 

3.0 R + JX 

17.0 

Osc. 

81.0 

Div. 

Div. 

3.B R + JX 

26.5 

Div. 

Div. 

Div. 

Div. 

4.0 R + JX 

76.0 

Div. 

Div. 

- Div. 

Div. 















TABLE - 3.3 


Effect of Large Line RaLtio on Convergence Ciype-b Simulation^ 


AD IEEE - 14 Buss System : 


111 -eondi tl oni ng 

'Cr A-S#* 

NurnhwAr of I i Arr all onsi 

R + JC1-/DX 

FDLP 

FIX-R 

NFC«-I 

NFDR-I M 

NFDR-I I 

R + jO. 75X 

4.B 

7. B 

0. 5 

0.5 

Div. 

R + J0.50X 

7.0 

8.0 

0.0 

8.5 

Div. 

R + JO. 2BX 

17. 0 

IS. B 

18.5 

0.0 

Div. 

R + JO. aox 

a3. B 

83.6 

ae. 0 

10.5 

Div. 

R + JO.lOX 

88. 0 

Osc. 

04.5 

1 

15.0 

Div. 

BD IEEE - 30 Bus 

System : 





111 -condi t 1 onl ng 

Number of Iterations 

R + jCl-^X 

FIX-P 

FCa-R 

NFDR-I 

NFDR-I M 

NFC«-II 

R + JO. 7SX 

I 5. 0 

7.0 

11.5 

11.5 

Div. 

R + JO. BOX 

7. 0 

O. 0 

0.5 

8.5 

Div. 

R + JO. SBX 

17. O 

10.0 

14.5 

8. 0 

Div. 

R + JO. aox 

84. B 

37. B 

ao.o 

0.0 

Div. 

R + JO. lOX 

Ofic. 

Obc. 

73.0 

14.0 

Div. 














3. 7 Discussions and Conclusions 


In this chaptor , the fast decoupled load flow methods in 
polar and in rectangular coordinates are presented. Three new 
formulations for FDLR are also suggested. All these methods 
have been tested on IEEE ~ 14 Bus and 30 Bus test systems both for 
the base case as well as ill-conditioned situations simulated for 
large R/OC ratio of lines. 

All the five methods CFDLP» FDLR» NFDR-I . NFDR-IM and 
NFI>R-II> utilize constant Jacobian's of similar size and thus have 
identical memory requirements. All these methods exhibit a 
geometric convergence pattern. 

Based on the results shown in Tables 3.1 to 3. 3» the relative 
computational performance of different methods are sumtmarized 
below : - 

1. From the CPU time tabulated in the present chapter and 
Chapter -2. it is clearly evident that all the five versions 
of fast decoupled load flow methods are faster than 
Newton-Raphson methods in polar as well as in rectangular 
coor di nates . 

2. Although, FDLR takes more number of iterations for 

convergence as compared to FEX-P, the total CPU time is 
comparable with FISLP model. 

3. The computational time required per iteration of FDLR is less 
than the corresponding Itefation time for FDLP model. 


4. Th« proposed NFDR-I mod*! though t*k»fi mor# number of 

lteraition« for b*se case» *.s compered to existing FEM-R, it 
hes better convergence chereeteri sties eapecielly for 

ill-conditioned ceses simuleted by inereesing the resistance 
of lin es. 

5. Thm proposed NFDR-II model takes maximum number of iterations 
for base ease and also diverges for ill-conditioned cases. 
Hence it is not recommended for load flow solution. 

6. The proposed modified version NFC^— IM model takes least 
number of iterations for ill -eondlton€Kl eases C simulated by 
decreasing the line reactances^ as compared to all other 
versions of fast decoupled load flow models. 



CHAPTER - IV 


CONCLUSION 


The work* carried out in ihi* theal* have been aimed at 
comparing different existing load flow methods, viz. 
Newton-Raphson method & fast decoupled load flow methods in polar 
and in rectangular co-ordinates as v*ell as developing improved 
solution technique* for power system load flow studies. Three new 
formulations of fast decoupled load flow methods In rectangular 
co-ordinates have been presented In this thesis. The inqportanee 
of any load flow solution depends largely up^on its merits 
regarding reliability, convergence characteristics, solution time 
and memory requirements. The above methods differ, mostly, in 
their memory and computation time requirements. Keeping this in 
view a comparative study of the aforementioned methods has been 
made. 


The potential of the existing NR and FDLF methods in polar 
and rectangular co-ordinates and new algorithms of FE>LF methods in 
rectangular co-ordinate* have been tested on IEEE-14 bus and 30 
bus systems both for well behaved and ill-conditioned cases. From 
the comparison of the test results of all these methods presented 
in Chapter 2 and 3, it is concluded that, 

CID Newton-Raphson load flow methods in polar and rectangular 
co-ordinates are found to 1^ highly reliable for solving well 
behaved and ill-conditioned systems. The number of 

iterations taken by these methods are almost Independent of 


system size and kind. 
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C 2 :> th* two N*wtori-R*phson method*, th* pol*r version 

CNRP> requires less memory then the recienaular versionCNRRS 
whereas the formulation in rectangular version is found to be 
slightly faster than it* ppolar counterpart. 

C33 Fast decoupled load flow method in polar co-ordinates and all 
its four version* in rectangular co-ordinates are found to be 
much faster than Newton -Raphson methods in polar as well as 
in rectangular co-ordinate* and also have lesser memory 
requirements as compared to the two NR methods. 

C45 For base case studies although the fast decoupled load flow 
™,thod. in r«=C.ngul.r co-ordin.t« CFDLR. t.Ke »one nunber 
of ItT.tlon. .« conp.r«i i*-- fonmul.Unn in pniT 

eo-ordin«t.s CFDUPS. ih» ''' eo«p»r»bl« with 

th. FDUP In f«t th. CPU ti». r«,uir«J pT it.r.ti»n 

by FDLR BOd.l i. l.«« <.h«" “"y ™''‘ ’ 

CE> Amongst th. thr.. n.v( formul.tlon. of f..*- dKOupl 
flow inod.1* in p.ct.ngul*r co-ordinAt« 

Ch.pt.r 3. NFDR-I .nd MFDR-IM vroion .r. cl«rly »op-rlor to 
th. .jd.ting FOLR mod«l «p«:i*lly solving ill condltl 

systsms. Amongst th. two. NFI»-X- 

oonv.rg.nc. p.tt.rn which is .y.n Fup.rioP to FDLP «od.l in 

certain situations. 

audios c.nni«i odt in this th.si. on f.*!- ^ 

flow „.thods ifi-«* - “ 

limits of generators etc. Iri ae 
considering the reactive power limitfS o a , * , 

..t a comparative study of relative 

the main emphasis was to present 

.«t.lna and new load flow models, 
performances of various existl g 


However, * pr*ciic*l load flow program should incorporate various 
syetem constraints such as generator reactive power limits, limits 
on transformer tappings, voltage limits at load buses etc. 
Efficient modelling of these constraints without hampering the 
constant Jacobian feature of the fast deoupled load flow methods 
has been an area of research till recerni past. Hence, the present 
work may be extended in future to incorporate this important 
feature in various fast decoupled load flow programs, which would 
further enhance its utility. 
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APPENDIX - A 


DATA FOR IEEE"-14. BUS TEST SYSTEM CAT 100 MVA BASD 

The IEEE-14 Bus system is shown in fig. A. The system data 
is taken from ref. C93 and buses renumbered. The data is compiled 

in the following manner i 

1) Bus data 

E) Line data 

3) Transformer data 

4) Shunt capacitor data 

5) Voltage controlled or P-V bus data. 

A-1 Bus Data 


Bus No. Generation Load 



H13 

fwar 

nw ” 

MVAR 

1 

0.0 

0.0 

0.0 

0.0 

E 

40.0 

0.0 

El. 7 

IE. 7 

3 

0.0 

0.0 

94. E 

19.0 

4 

0.0 

0.0 

11.2 

7.5 

5 

0.0 

0.0 

0.0 

0.0 

6 

0.0 

0.0 

47.8 

-3.9 

7 

0.0 

. 0.0 

0.0 

0.0 

8 

0.0 

0.0 

7.6 

1 .6 

9 

0.0 

0.0 

29.5 

16 . 6 

10 

0.0 

0.0 

9.0 

5.8 

11 

0.0 

0.0 

3.5 

1*8 

IE 

0.0 

0.0 

6.1 

1.6 

13 

0.0 

*0.0 

13.5 

5 • 8 

14 

0.0 

0.0 

14.9 

5.0 



n.-i! Line Data 


Duse f. 

Line ^4o. i^root' To 

l ine l»f(p»»danc<? 

R in p.u. X in p.u. 

Half Line 
Charging 

susceptance 

(p.u) 

1 

1 

2 

0.01938 

0.03917 

0.0264 

2 

1 

8 

0.05403 

0.22304 

0.0246 

3 

2 

3 

0.04699 

0.19797 

0.0219 

4 

2 

6 

0.05811 

0.17632 

0.0187 

5 

2 

8 

0.05695 

0.17388 

0.0700 

6 

3 

6 

0.06701 

0.17103 

0.0173 

7 

4 

11 

0.09498 

0.19890 

0.0 

S 

4 

12 

0.12291 

0.25581 

0.0 

9 

4 

13 

0.06615 

0.13027 

0.0 

'10 

6 

7 

0.0 

0.20912 

0.0 

11 

6 

8 

0.01335 

0.04211 

0.0064 

12 

6 

9 

0.0 

0.55618 

0.0 

13 

7 

5 

0.0 

0. 17615 

0.0 

14 

7 

9 

0.0 

0.11001 

0.0 

15 

8 

4 

0.0 

0.25202 

0.0 

16 

9 

10 

0.03181 

0.08450 

0.0 

17 

9 

14 

0.12711 

0.27038 

0.0 

18 

10 

11 

0.08205 

0. 19207 

0.0 

19 

12 

13 

0.22092 

0. 19988 

0.0 

20 

13 

14 

' 0.17093 

0.34802 

0.0 



A- 

-3 Transformer Data 

Transformer 

From 

Buses setting 

To 

1 

6 

7 0.978 

2 

6 

9 0.969 

3 

8 

4 0.932 


A--4 

t 

Shunt Capacitor Data 



Bus No, 


Susceptance p.u. 

9 


0.190 


I 



A-5 Voltage-Controlled Bui Dsta (P-V buses) 




Bus No. 

Vo 1 tage 
magnitude p. u. 

Reactive 
Minimum NVAR 

Power Limits 

Maximum mVaR 

2 

1.045 

1 

O 

• 

O 

50.0 

3 

1-0 10 

c 

m 

o 

0.0 

4 

1 .070 

0 

m 

1 

34.0 

5 

1,090 

- 6.0 

o 

m. 


Slack bus voltage = 1,06/0^ 



AF» Pr.NDiX ~ B 

DATA FOR CEE- 30 BUa TEBT 9V8TEM (AT 100 MVA BA3E) 


S -1 Ehus Data 


Bus No . 

G «» n»»ration 

{y(|3 

HW 

Load 

MVAR 

1 

0.0 

* 0.0 

0.0 

0.0 

2 

40.0 

0.0 

21.7 

12.7 

3 

0.0 

0,0 

30.0 

30.0 

4 

0-0 

0.0 

0 . 0 

0.0 

5 

0.0 

0.0 

94.2 

19.0 

6 

0,0 

0.0 

0.0 

0.0 

7 

0.0 

0.0 

22.8 

10.9 

8 

0.0 

0.0 

2.4 

1.2 

9 

0.0 

0.0 

0.0 

0.0 

10 

0.0 

0,0 

5.8 

2.0 

11 

0.0 

0.0 

7.6 

1.6 

12 

0.0 ’ 

0.0 

11.2 

7.5 

13 

0.0 ' 

0,0 

0.0 

0.0 

14 

: 0.0 

0,0 

6.2 

1.6 

15 

0.0 

0.0 

8.2 

2.5 

16 

0.0 

0.0 

3.5 

1.8 

17 

0.0 

0.0 

9.0 

5 .8 

18 

0-0 

0.0 

3.2 

0.9 

19 

0.0 

0.0 

9.5 

3.4 

20 

0.0 

0.0 

2.2 

0.7 

21 

0.0 

0.0 

17.5 

11.2 

22 

0.0 

0,0 

0.0 

0.0 

23 

0.0 

0.0 

3.2 

1*6 

24 

0.0 

0.0 

8.7 

6 • 7 

25 

0.0 

0.0 

0.0 

0.0 

26 

0,0 

0.0 

. 3.5 

2.3 

27 

0.0 

0.0 

0.0 

0.0 

28 

0.0 

0.0 

0.0 

0.0 

29 

0.0 

0.0 

2.4 

0.9 

30 

0.0 

* 0.0 

10.6 

1.9 


( I'l# tmtf\ 


Line No, 


Fuse*! 

Frofir“To 


a 

a 

11 

11 

5 

13 

13 
7 
7 

3 
9 

10 

4 

10 

12 

h 

14 

15 
. 16 
' 15 

17 

18 
,19 
20 
20 
17 
21 
22 
22 

23 

24 

24 

25 

26 
27 

27 

29 

30 
30 

28 
28 


0.0192 

0.0452 

0.0570 

0.0132 

0.0472 

0.0581 

0.0119 

0.0460 

0.0267 

0.0120 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

• 0.1231 
, 0.0662 
0.0945 
0.2210 
0.0824 
0.1070 
0.0639 
0.0340 
0.0936 
0.0324 
0.0348 
0.0727 
0.0116 
0.1000 
0.1150 
0.1320 
0.1885 
0.2344 
0,1093 
0.0 
0.2198 
0.3202 
0.2399 
O . 0636 
0.0169 


X ' ifi ' p . u . 

Half Line 
Charging 

susceptani 

( p . u ) 

0.0375 

0.0264 

0. 1852 

0.0204 

0.1737 

0.0184 

0 , 0379 

0.0042 

0. 1983 

0.0209 

0. 1763 

0.0187 

0.0414 

0.0043 

0.1160 

0.0102 

0.0820 

0.0085 

0.0420 

0.0045 

0.2080 

0.0 

0.5560 

0.0 

0.2080 

0.0 

0. 1100 

0.0 

0.2560 

0.0 

0.1400 

0.0 

0.2559 

0.0 

0.1304 

0.0 

0.1987 

0.0 

0.1997 

0.0 

0.1923 

0.0 

0.2185 

0.0 

0. 1292 

0.0 

0.0680 

0.0 

0.2090 

0.0 

0.0845 

0.0 

0.0749 

0.0 

0.1499 

0.0 

0.0236 

0.0 

0.2020 

0.0 

0. 1790 

0.0 

0.2700 

0.0 

0.3292 

0.0 

0.3800 

0.0 

0.2087 

0.0 

0.3960 

0.0 

0.4133 

0.0 

0.6027 

0.0 

0.4533 

0.0 

0.2000 ■ 

0.0214 

0.0599 

0.0065 


P-y ,U' «, f ‘I f m»> r 


Transformer 


U y t* ^ 

Tap fatting 

FrofT 

'"'To 

1 

11 

12 

0.932 

Z 

13 

f * 

0.978 

, 3 

13 

10 

0.989 

4 

28 

27 

0.968 


B~4 

Shynt Capacitor Data 



Bus No. 

Susceptance p.u. 

10 

0.190 

24 

0.043 

B~3 Voltage-Controlled Bus Data CP~V buses) 

Bus No. 

Voltaoe Reactive Power Limits^ 

magnitude p.u. Minimum MVAR naKimum MVAR 


2 

1.045 

- 40.0 

30.0 

3 

1.010 

- 10.0 

40.0 


1.0B2 

- 6.0 

24.0 ; 

5 

1.010 

. 40.0 


6 

1.073 

- 6.0 

24.0 



Slack bus 

vol taiQ©, 1 n0&/0 ■■ ;j 


i 
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